Report  No.  91015 


AD-A241  596 


UNLIMITED 


Ami  J  Ks 


Report  No.  91015 

ROYAL  SIGNALS  AND  RADAR  ESTABLISHMENT, 

MALVERN 


C 


M.  .1 

ELECTE 

OCT  1  5  1991 


A  REVIEW  OF  POSITIVE  ION  SENSITIVITIES 
FOR  THE  SIMS  ANALYSIS  OF  CMT 


Author:  R  Holland  &  GW  Blackmore 


This  document  has  boon  approved 
for  public  release  and  sale;  its 
distribution  is  unlimited. 


iiiiiiiif 


PROCUREMENT  EXECUTIVE,  MINISTRY  OF  DEFENCE 

RSRE 

Malvern,  Worcestershire. 


May  1991 


UNLIMITED 

91  loll  o?7 


01067^1 


CONDITIONS  OF  RELEASE 


303852 


DRIC  U 


COPYRIGHT  (c) 
1988 

CONTROLLER 
HMSO  LONDON 


. . .  DRIC  Y 

Reports  quoted  are  not  necessarily  available  to  members  of  the  public  or  to  commercial 
organisations. 


DEFENCE  RESEARCH  AGENCY  ELECTRONICS  DIVISION 


Report  91015 

TITLE:  A  REVIEW  OF  POSITIVE  ION  SENSITIVITIES  FOR  THE  SIMS 

ANALYSIS  OF  CMT 

AUTHORS:  R  Holland  (Consultant)  and  G  W  Blackmore 
DATE:  May  1991 


SUMMARY 

Measurement  of  the  relat.ve  intensities  of  positive  ion  emission  from  CMT  under  SIMS 
analysis  provides  calibration  factors  f^r  the  high  sensitivity  analysis  of  impurities.  The 
Relative  Sensitivity  Factor  (RSF)  values  found  by  RSRE,  mainly  from  uniformly  doped 
samples,  compared  well  with  those  obtained  late:  by  other  workirs  using  implants. 
Theoretical  models  for  sputtered  ion  emission  are  reviewed  in  detail  and  their 
applicability  to  the  CMT  data  is  considered,  particularly  in  relation  to  plots  of  log  (RSF) 
against  ionisation  potential.  It  is  suggested  that  more  than  one  mechanism  is  likely  to 
be  operative.  Comparison  with  the  results  for  other  semiconductors  suggests  that  the 
observed  spread  from  a  linear  relationship  might  be  attributable  to  a  low  implanted 
oxygen  concentration  resulting  from  a  high  sputter  rate. 
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INTRODUCTION 


The  use  of  Secondary  Ion  Mass  Spectrometry  (SIMS)  analysis  has  been  an  essentia!  factor 
in  the  progress  of  the  RSRE  programme  to  develop  infrared  sensitive  materials  and 
devices  based  on  high  purity  cadmium  mercury  telluride  (CMT).  SIMS  identifies  and 
counts  the  ionised  particles  sputtered  from  a  surface  bombarded  by  a  focussed  beam  of 
ions  of  several  keV  energy.  The  resulting  analysis  can  cover  all  elements  with 
unrivalled  sensitivity  (down  to  1013atoms/cm°  in  favourable  cases)  and  can  be  highly 
localised,  typically  to  less  than  a  micron  in  width  and/or  depth.  It  is  thus  particularly 
suitable  for  the  identification  and  location  of  dopants  and  impurities  in  epitaxial  layers 
and  device  structures. 

While  the  elemental  composition  of  the  sputtered  material  is  representative  of  that  of 
the  parent  material,  SIMS  can  only  detect  the  ionised  components.  The  degree  of 
ionisation  varies  widely  between  elements  -  eg  between  alkali  metals  and  gold  -  by  as 
much  as  5  decades  and  is  also  dependent  on  the  bulk  chemistry  of  the  parent  surface, 
such  as  the  matrix  composition  and  its  state  of  oxidation.  Attempts  to  increase  the 
degree  of  ionisation  by  the  use  of  lasers  or  plasma  discharge  have  had  only  limited 
success.  Fortunately  experience  shows  that  under  coasrant  experimental  conditions  the 
degree  of  ionisation  remains  essentially  constant  over  wide  ranges  of  concentration, 
from  parts  per  billion  to  a  few  percent,  so  that  variations  in  the  number  of  ions 
detected  accurately  reflect  proportionate  changes  of  composition  with  depth.  In  order 
to  quantify  such  data,  however,  the  analyst  requires  calibration  factors  for  each 
element  of  interest  in  each  matrix. 

Singly  charged  monatomic  and  diatomic  ions  predominate  in  SIMS  spectra.  The  degree 
of  ionisation  of  most  electropositive  elements  is  greatly  enhanced  by  surface  oxidation, 
either  by  use  of  02*  or  O'  primary  ions  or  by  back-flooding  the  sample  chamber  with 
oxygen  at  pressures  up  to  lO'^orr.  Similarly  the  production  of  negative  ions  of 
electronegative  elements  is  enhanced  by  caesium.  The  present  report  is  concerned  with 
the  establishment  and  correlation  of  reliable  calibrations  for  singly  charged  positive 
monatomic  ions  sputtered  from  CMT  by  oxygen  primary  ions.  Particular  attention  is 
paid  to  metallic  and  metalloid  elements  with  ionisation  potentials  in  the  range  5-10eV. 
For  elements  with  higher  ionisation  potentials  caesium  primary  ions  are  usually  used. 

During  the  past  decade  calibration  factors  have  been  measured  in-house  using  doped 
crystals,  ion  implants  and  diffused  radio-tracers.  These  can  now  be  compared  with  data 
published  by  other  workers  and  the  predictions  of  theoretical  models  in  order  to 
establish  the  degree  of  confidence  in  their  use,  to  tdentify  those  values  which  need 
further  examination  and  to  predict  likely  values  for  elements  which  may  not  have  been 
measured  previously. 

The  simplest  correlation  betwseen  positive  secondary  ion  sensitivities  in  a  given  matrix 
is  offered  by  the  well-known  inverse  relationship  between  ionisation  potential  and  the 
logarithm  of  the  sensitivity.  The  theoretical  basis  for  this  relationship  and  the  other 
factors  which  have  been  introduced  to  improve  it  are  discussed  later.  Lapides  and  co¬ 
workers1  ’2  published  sensitivity  factors  for  10  positive  and  10  negative  secondary  ions 
from  CMT.  The  values  for  8  of  the  positive  ions  showed  a  reasonable  correlation  with 
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ionisation  potentials  while  those  for  2  elements,  Mg  and  Be,  departed  from  the  trend 
line  by  a  decade.  Wilscn3,4  published  values  for  42  elements  and  showed  a  similar 
correlation,  Mg  and  Be  again  lying  off  the  line. 

As  the  observed  ion  yieiu  is  dependent  on  instrumental  factors  which  may  vary  between 
measurements,  a  suitable  matrix  isotope  is  usually  chosen  as  an  internal  standard.  The 
term  'Relative  Sensitivity  Factor'  (RSF),  now  commonly  used5,  refers  to  the 
concentration  of  the  impurity  or  dopant  which  would  (assuming  linear  response)  give  the 
same  number  of  secondary  ions  as  the  matrix  element.  Measurements  based  on 
particular  isotopes  should  therefore  be  divided  by  the  appropriate  isotope  fraction  to 
give  the  total  yield  for  the  measured  or  matrix  element.  The  term  'Relative  Sensitivity 
Factor'  is  a  misnomer  as  it  varies  inversely  with  the  true  sensitivity  of  the  element. 
Nevertheless,  it  will  be  used  here  as  it  is  coming  into  common  usage5.  The 
concentration  c  is  given  by: 

c  =  RSF.(yiN)  .  (fN/y 

where  I*  and  IN  aie  the  secondary  ion  currents  for  isotopes  M  and  N  of  the  element  and 
matrix  and  fM  and  fN  are  the  respective  isotope  ratios. 

EXPERIMENTAL 

SIMS  measurements  were  made  using  the  Cameca  IMS-3f  ion  microscope  formerly  used 
at  AMTE,  Holton  Heath,  and  now  installed  at  RSRE.  A  beam  of  8keV  02+  primary  ions 
was  rastered  over  a  500/im  square  area  of  the  sample.  The  secondary  ion  optics  were 
adjusted  to  measure  the  emission  from  a  central  circular  area  of  60^m  diameter.  The 
125Te  isotope  was  chosen  as  a  reference  standard  as  it  is  free  from  molecular  ion 
interferences. 

Doubly  charged  ions  of  Cd  and  Te  were  found  to  interfere  extensively  with  the 
measurement  of  singly  charged  ions  with  the  same  mass/charge  ratio.  This  affected  the 
transition  elements  Mn,  Fe,  Ni,  Cu  and  Zn.  Because  of  the  low  energy  of  the  interfering 
ions,  it  was  shown  that  this  effect  could  be  reduced  by  displacing  the  energy  band-pass 
by  lOeV6,7.  The  origin  of  this  energy  scBle  was  arbitrarily  taken  from  the  energy  at 
which  the  intensity  of  AF  ions  used  to  set  up  the  instrument  was  a  maximum.  All 
measurements  were  made  with  a  band-pass  of  lOeV  centred  on  this  displacement.  The 
energy  distributions  of  most  of  the  elements  of  interest  were  found,  in  general,  to 
parallel  that  for  the  Te+  ions  used  as  a  matrix  reference  (Fig  1).  Consequently  the 
choice  of  band-pass  made  little  difference  to  the  measured  RSFs,  apart  from  the 
removal  of  interferences.  One  exception  is  As,  as  discussed  later. 

Several  doped  crystals  were  prepared,  the  most  satisfactory  one  being  made  by  melting 
powdered  CMT  intimately  mixed  with  finely  powdered  tellurides  (specially  made  at 
BDH,  Poole)  of  the  range  of  elements  of  Interest.  The  particles  all  passed  through 
30^an  mesh  nylon  selves8.  Where  tellurides  were  unsuitable,  appropriate  alternative 
stable  compounds  such  as  cadmium  salts  were  incorporated.  The  resulting  crystal 
proved  to  be  of  remarkably  uniform  composition,  whereas  similarly  prepared  crystals 
of  cadmium  telluride  were  less  satisfactory  because  of  segregation  of  some  dopants  - 
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possibly  because  of  lower  solubilities.  The  dopant  concentrations  were  measured  by 
flameless  atomic  absorption  and  by  ICP  using  standard  solutions  for  calibration. 
Adjacent  slices  of  the  crystal  were  used  for  these  analyses  and  for  SIMS  measurements. 
Not  all  dopants  were  measured  but  the  agreement  of  the  analyses  with  the  amount  of 
dopants  used  was  so  consistant  that  the  few  elements  not  measured  were  assumed  to  be 
present  at  the  concentration  used  in  crystal  preparation. 

Ion  implants  were  also  used  to  establish  RSF  values  using  the  technique  of  Leta  and 
Morrison9  in  which  the  integrated  SIMS  count,  after  subtracting  the  background,  is 
compared  with  the  total  number  of  ions  implanted.  Applications  of  this  technique  were 
limited  by  the  restricted  availability  of  suitable  implants.  The  range  of  implanted  ions 
must  be  sufficient  for  them  to  be  clearly  distinguished  from  surface  contamination, 
discrimination  being  especially  difficult  for  fast  diffusing  ions10.  Doubts  were  also 
entertained  as  to  whether  damage  caused  by  the  implantation  wouid  affect  the  results. 

For  Au  and  As,  additional  data  was  obtained  by  collaboration  with  studies  at  the 
University  of  Southampton  in  which  profiles  of  the  diffusion  of  radioactive  tracers  into 
CMT  were  measured  using  an  etch  technique11'12.  Comparison  of  these  profiles  with 
those  obtained  by  SIMS  gave  revised  values  for  the  RSFs. 

RESULTS 

Between  four  and  six  measurements  were  made  for  each  element  and  were  generally 
found  to  be  within  a  factor  of  2  of  the  mean,  readings  markedly  dissimilar  to  the  mean 
being  rejected.  Where  implants  were  available,  the  results  were  generally  found  to  be 
in  the  same  range.  As  the  values  for  different  elements  ranged  over  5  decades,  an  error 
factor  of  2  was  regarded  as  acceptable.  The  mean  values  are  listed  in  Table  1,  together 
with  those  published  later  by  Lapides  et  al  and  by  Wilson.  Fig  2  shows  a  plot  of  log 
(RSF)  versus  ionisation  potential  (in  eV)  and  Fig  8  is  the  corresponding  plot  for  Wilson's 
data.  Of  the  19  common  elements,  16  showed  differences  on  a  log  scale  of  0.3  or  less 
(ie  a  factor  of  2)  The  points  for  the  3  elements  which  differed  considerably  are  marked 
with  circles  in  Figs  2  and  8  and  the  two  values  obtained  on  diffusion  specimens  are 
marked  with  crosses  (x)  and  their  values  in  Table  1  are  enclosed  in  brackets. 

DISCUSSION 

Most  of  RSRE's  values  for  the  RSFs  listed  in  Table  1  were  obtained  in  1980  and  the 
degree  of  agreement  with  those  reported  by  other  workers  in  1984  and  1988  is 
remarkable  in  view  of  the  very  different  methods  employed.  The  results  reported  by 
Lapides  et  al  and  by  Wilson  were  all  obtained  from  Implants,  whereas  RSRE's  results 
were  nearly  all  based  on  doped  crystals.  Furthermore,  it  is  most  unlikely  that  the 
former  workers  used  the  somewhat  arbitrary  choice  of  energy  band-pass  used  by  RSRE. 
Their  measurements  were  made  by  Evans  Associates  and  the  experimental  conditions 
were  not  reported.  Evans'  usual  practice  had  been  to  use  a  wide  band-pass  of  some 
120eV  with  an  offset  of  some  30eV  when  mass  interferences  are  suspected.  The 
agreement  is  a  fortunate  consequence  of  the  similarity  between  the  energy  distributions 
of  the  sputtered  ions  and  that  of  Te*  used  as  a  reference  element.  It  suggests  that 
these  RSFs  may  be  used  with  confidence  and  could  also  be  used  as  a  basis  for  any 
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generalised  approach  to  establish  a  correlation  between  different  elements.  It  also 
suggests  that  damage  caused  by  ion  implantation  does  not  affect  SIMS  calibrations  in 
this  matrix. 

Tne  3  elements  for  which  Wilson's  and  RSRE's  values  differ  widely  are  As,  Sb  and  Cu. 
The  number  of  measurements  on  which  Wilson's  values  in  Table  1  were  based  averaged 
5  per  element.  Those  for  As+  and  Sb+  were  based  on  one  each,  although  B+  was  based 
on  16  readings  and  In*  on  13.  It  can  be  seen  from  Fig  1  that  the  energy  distribution  for 
As*  differs  in  shape  from  luat  of  the  reference  Te*,  the  maximum  of  the  population 
curve  being  displaced  to  lower  energy  by  some  lOeV  and  the  fall-off  towards  higher 
energies  is  steeper.  Consequently,  in  view  of  the  different  energy  windows  used, 
agreement  between  the  2  laboratories  could  not  be  expected.  The  value  based  on 
diffusion  profiles,  involving  many  measurements  using  the  maximum  energy  band-pass 
of  the  instrument,  are  probably  the  most  reliable.  The  concentration  of  Sb  in  the  RSRE 
doped  sample  was  not  measured  directly  but  based  on  the  amount  added.  The  SIMS 
count  was  very  low  and  the  resulting  RSF  must  be  suspect.  Doubt  may  also  be  cast  on 
RSRE's  value  for  Cu*  as  the  interference  from  Te2*  at  mass  63  and  65  is  only  partially 
reduced  by  the  voltage  offset  used13.  Nor  can  one  choose  between  Wilson's  value  and 
that  of  Lapides,  which  differed  by  a  factor  of  5.6.  It  was  hoped  that  comparison  with 
the  RSFs  of  other  elements  using  semi-empirical  correlations,  preferably  supported  by 
theoretical  models  of  the  ionisation  mechanism,  would  help  to  resolve  these  anomalies. 

There  have  been  many  attempts  to  formulate  a  general  quantitative  model  for 
secondary  ionisation.  The  most  successful  in  terms  of  matching  known  data  and 
predicting  new  values  in  a  variety  of  matrices  was  the  so-called  Local  Thermal 
Equilibrium  (LTE)  model  proposed  by  Andersen  and  Hinthorne14.  They  suggested  that 
the  sputtering  site  could  be  likened  to  a  plasma  in  which  atoms,  ions,  electrons  and 
molecular  ions  were  in  local  equilibrium.  They  therefore  used  the  Saha-Eggert  equation, 
which  successfully  predicts  ion  intensities  in  flame  and  plasma  spectra,  to  predict 
secondary  ion  concentrations.  This  takes  the  form: 

Y  =  A.(BVB0)  .  exp(-E/kT) 

where  Y  is  the  ion  yield  for  a  particular  element,  A  is  a  constant  for  the  matrix,  B°  and 
B1  are  the  electronic  partition  functions  of  the  unionised  and  ionised  atoms  of  the 
element  at  temperature  T,  k  is  Boltzman's  constant  and  E  is  the  ionisation  potential. 
A  similar  relation,  with  electron  affinity  replacing  ionisation  potential,  applied  to 
negative  ions.  These  relations,  extended  to  include  oxide  ions,  formed  the  basis  of  a 
commercial  computer  programme  called  CARISMA  for  use  with  the  ARL  ion 
microprobe.  Inter-laboratory  exercises  organised  by  NRL  using  standardised  glasses  and 
steels15  showed  considerable  agreement  usually  within  a  factor  of  2. 

Unfortunately  there  was  no  physical  basis  for  this  model,  despite  its  success.  The 
'temperature'  T  was  unrealistically  high,  usually  in  the  range  3000  to  10,000  degrees, 
and  differed  for  positive  and  negative  ion  emissions.  Morgan  and  Werner16,17,  using  an 
empirical  approach  and  accepting  these  limitations,  successfully  used  the  relationship 
in  simplified  form,  plotting  log  (Y.B°/B1)  against  E  to  find  the  best  straight  line.  The 
3lope  gives  the  term  1;T,  but  the  partition  functions  B  are  a  function  of  T  and  an 
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iterative  procedure  is  needed  to  find  the  best  fit.  B  is  defined  as  the  sum: 

Y  exp(-e/kT) 

over  all  electronic  energy  levels  e1  and  gives  a  weighting  for  rh<=>  inclusion  of 
electronically  excited  states.  Values  for  B  at  temperature  T  can  be  calculated  from 
published  data18.  Bearing  in  mind  that  RSF  is  reciprocally  related  to  Y,  Fig  3  shows  a 
plot  of  log  [(RSFKbVb0)  using  the  mean  log  (RSF)  values  of  Table  1.  The  slope  gives 
5500°K  for  the  value  of  T.  Values  of  bVb°  at  this  temperature  are  listed  in  Table  2. 
Fig  3  shows  only  a  slight  improvement  in  linearity  over  Fig  1,  the  shift  of  Mg  into  the 
main  family  being  offset  by  the  loss  of  linearity  among  elements  of  lower  ionisation 
potential  and  displacement  of  the  point  for  Ca. 

Morgan  and  Werner  introduced  the  square  root  of  the  atomic  weight  thus: 

Y  =  A.m'4  .  (BVb°)  .  exp(-E/kT) 

explaining  that  it  was  a  composite  factor  taking  into  account  the  mass  discrimination 
of  the  instrument  transmission  and  detection  as  well  as  in  ion  formation,  but  mainly 
because  it  improved  the  fit.  The  fit  was  also  improved  by  using  the  sum  of  ion  and 
oxide  ion  yields.  Use  of  the  term  m~4  did  not  improve  the  fit  of  the  present  data. 
Slodzian19,  studying  isotopic  discrimination  in  SIMS,  generalised  the  term  to  m"Q,  with 
o  lying  usually  in  the  range  0  5  to  2.  There  was  no  besr  value  of  a  to  fit  the  present 
data  and  its  use  appears  artificial. 

The  measurement  by  Rudat  and  Morrison20  of  correction  factors  to  convert  the  signals 
from  photo-multipliers  into  true  ion  currents  is  c:  more  practical  significance.  The 
factors  vary  from  element  to  element  but  are  related  more  to  the  chemistry  of  the 
element  than  to  its  atomic  mass.  The  data  listed  in  Table  2  shows  by  what  factor  the 
true  current  is  enhanced.  The  variation  between  these  factors,  however,  is  not 
sufficient  to  reduce  the  scatter  of  RSFs  significantly  (Fig  4). 

Despite  its  man>  successes  in  predicting  secondary  ion  sensitivities  in  a  variety  of 
matrices,  Andersen  and  Hinthorne’s  LTE  model  is  inconsistent  with  current  ideas  of  the 
mechanism  of  sputtering.  The  incoming  02*  ion  is  believed  to  initiate  a  collision 
cascade  while  penetrating  several  atomic  layers  and  coming  to  rest  as  implanted  atoms. 
Computer  simulations  show  that  the  cascade,  which  dissipates  in  about  10~14  second, 
causes  some  atoms  to  be  ejected  from  the  surface  layer  by  impacts  from  within  the 
solid21.  This  timescale  is  insufficient  by  several  orders  of  magnitude  for  thermal 
equilibrium  to  be  established.  A  primary  current  of  1/iA  corresponds  to  6xl012  impacts 
per  second  so  that  even  if  they  all  occured  at  the  same  spot  they  would  still  be  isolated 
events. 

The  surface  region  of  the  specimen  must  contain  a  high  proportion  of  implanted  oxygen 
which  is  progressively  sputtered  as  it  is  exposed  on  the  surface.  To  maintain 
equilibrium  each  sputtering  event  should,  on  average,  involve  the  loss  of  2  oxygen 
atoms.  With  a  primary  ion  current  of  1/iA  on  l/4mm2,  the  rate  of  erosion  of  CMT  found 
in  the  present  work  is  5/im/hour,  which  ccrrespcr.dc  to  1.6  matrix  atoms  pci  impact,  in 
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addition  to  the  2  oxygen  atoms.  Changing  the  primary  ion  current  would  not  affect  this 
proportion  provided  the  erosion  rate  remains  proportional  to  the  current.  It  also  means 
that  the  oxygen  concentration  in  the  surface  layers  is  1.25  atoms  per  matrix  atom. 
These  oxygen  atoms  might  well  occupy  the  vacancies  caused  by  the  inevitable 
preferential  loss  of  mercury  as  well  as  interstitial  sites.  Alternatively,  allowing  for  the 
loss  of  mercury,  there  would  be  sufficient  oxygen  to  convert  all  the  remaining  matrix 
atoms  to  oxides,  Te02  and  CdO. 

Any  general  theory  of  the  ionisation  of  sputtered  particles  must  account  for  the  broad 
applicability  of  the  Saha-Eggert  relation,  the  variation  between  matrices  and  the 
enhancement  of  ion  yields  by  oxygen.  Morgan  and  Werner16  pointed  out  that  while  the 
energy  imparted  to  any  one  atom  is  random,  the  statistical  effect  should  be  the  same 
for  all  elements  provided  the  differences  in  cross-section  can  be  discounted.  The  energy 
will  be  used  partly  to  overcome  the  binding  forces,  partly  as  kinetic  energy  and  partly 
as  excitation  energy  distributed  among  the  quantum  levels.  The  probability  of  atoms 
acquiring  excitation  energy  ei  can  be  shown  by  statistical  mechanics  (see  Appendix  A) 
to  take  the  Boltzman  form  A  exp  (-fie,),  where  P  is  analagous  to  VkT  and  kT  is  a 
measure  of  the  mean  energy  available.  From  this  it  follows  that  the  ratio  of  ionised  to 
neutral  atoms  is 


BVB0  exp  (-0E) 

This  approximates  to  the  Saha-Eggert  equation  when  the  proportion  of  neutrals  exceeds 
that  of  the  ions. 

This  approach  is  grossly  oversimplified.  Apart  from  taking  no  account  of  tne  different 
binding  energies,  it  appears  to  over-emphasise  the  role  of  neutral  atoms  rather  than 
oxidec  and  oxide  ions.  It  is  probably  more  applicable  to  metallic  matrices  than  ionic 
lattices  and  is  possibly  more  valid  for  the  metal  impurities  than  for  alkali  metals  or 
alkaline  earths.  One  would  not  expect  to  find  neutral  alkali  metal  atoms  in  such  an 
oxidising  environment.  In  recent  years  it  has  become  customary  to  omit  the  partition 
functions  and  to  plot  log  (RSF)  directly  against  ionisation  potential  on  an  empirical 
basis22. 

Strictly,  as  ionisation  involves  donating  an  electron  to  the  surface,  E  should  be  replaced 
by  (E-0),  where  0  is  the  work  function.  This  would  be  included  in  the  constant  A. 
Deline  et  al23  proposed  that  the  matrix  effect  could  be  entirely  accounted  for  by  the 
effect  of  the  work  function.  It  is  increased  by  oxygen  coverage  thus  enhancing 
ionisation.  Under  oxygen  bombardment  the  oxygen  concentration  in  the  surface  layer 
will  be  greater  the  lower  the  sputter  rate  as  discussed  above.  However,  it  is  unlikely 
that  the  effect  would  be  the  same  for  all  elements,  eg  one  would  expect  a  greater 
enhancement  for  iron  and  nickel  than  for  gold. 

The  model  outlined  above  is  useful  but  is  unlikely  to  reflect  the  true  situation. 
Slodzian24  suggested  that  oxygen  chemically  reacts  with  the  impurity  elements  so  that 
sputtering  involves  breaking  oxide  bonds.  This  'bond-breaking'  model  is  now  widely 
accepted  though  it  is  difficult  to  quantify.  Lang  and  Yu25,26,27  have  tried  to  estimate 
the  probability  that  ions  leaving  the  surface  would  be  neutralised  while  traversing  the 
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surface  field  be.  ,eir  calculations  of  the  overlap  of  electron  orbitals  are  not  readily 
translated  in.o  practical  RSFs. 

It  is  evident  that  ionisation  is  a  complex  process  and  several  mechanisms  may  be 
operative  not  only  in  different  matrices  but  possibly  for  different  elements  in  the  same 
matrix.  In  ionic  solids  some  ions  may  be  ejected  directly  from  the  lattice.  In  other 
cases  neutral  atoms  may  be  ejected  in  highly  excited  electronic  states,  possibly  by 
promotion  of  electrons  from  inner  orbitals.  After  a  short  lifetime  an  electron  would 
then  be  ejected  by  an  Auger  process24  outside  the  surface  and  the  process  would  be 
independent  of  the  work  function.  Alternatively,  for  elements  which  readily  form 
oxides,  short-lived  oxide  molecules  with  high  degrees  of  rotational,  vibrational  and 
electronic  excitation  may  also  decompose  outside  the  surface: 

MO  -  M*  +  O' 

Examination  of  Fig  2  supports  this  idea  that  several  mechanisms  are  operative.  The 
spread  covers  two  decades  with  the  nobler  elements  Ag,  Pt,  Au  and  the  non-metals  B 
and  Si  forming  the  upper  boundary  of  RSFs  and  the  elements  Mg  and  Be  forming  the 
lower  boundary.  This  suggests  that  the  latter  elements  are  ionised  by  a  different 
process  than  the  former.  Now  that  tables  of  RSFs  are  available  for  other 
semiconductors,  it  is  interesting  to  compare  them  with  those  for  CMT.  Figs  5  to  8  show 
the  plots  for  Si,  GaAs,  InP  and  CMT.  In  each  case  there  is  a  general  trend  relating  log 
RSF  with  E,  the  spread  increasing  towards  higher  values  of  E  and  with  Mg  and  Be 
forming  the  lower  boundary.  Indeed  Mg  and  Be  appear  to  form  a  class  by  themselves. 

It  is  tempting  to  suggest  that  the  oxide  decomposition  mechanism  applies  to  Mg  and  Be. 
Whereas  an  excited  atom  needs  energy  E  to  ionise,  an  oxide  molecule  or  radical  needs 
D  ♦  E  -  Eo  where  D  is  the  heat  of  dissociation  and  E0  is  the  electron  affinity  of  oxygen 
(1.46eV).  If  D  is  already  acquired  during  the  sputtering  event  (eg  by  bonding  with  a 
moving  oxygen  atom),  the  effective  ionisation  energy  E  is  reduced  by  1.46eV,  which 
could  account  for  the  displacement  of  the  points  for  Mg  and  Be  on  the  graph  from  the 
main  trend  line.  The  high  concentration  of  implanted  oxygen  would  make  oxygen  atoms 
the  most  common  particles  in  the  cascade  and  their  low  mass  would  ensure  they  were 
the  most  mobile.  It  is  possible  that  this  type  of  mechanism  could  also  apply  even  if  the 
oxide  molecule  did  not  escape  the  surface  -  it  could  be  formed  transiently  by  oxygen 
atom  impact,  the  O  ion  remaining  in  the  surface.  There  is  no  obvious  way  of  checking 
this  theory. 

If  one  discounts  Mg  and  Be,  it  is  evident  that  the  spread  cf  the  plot  for  Ir.P  is  much 
greater  than  that  for  Si  and  GaAs  and  more  like  that  for  CMT.  This  is  illustrated  by 
the  spread  of  the  points  for  Mn,  Ni,  Sn  and  Ag,  all  with  similar  ionisation  potentials. 
The  differences  in  log  RSF  between  Mn  and  Ag  are  0.75  and  0.72  decades  in  S>  and  GaAs 
but  1.17  and  1.95  decades  in  CMT  and  InP  respectively.  Furthermore,  the  vertical 
spread  of  elements  tends  to  reflect  their  affinity  for  oxygen.  The  concentration  of 
oxygen  in  the  surface  layers  might  well  play  an  important  role  in  determining  this 
spread.  Morgan  and  Werner  showed  that  impurities  in  steels  only  fitted  their  modified 
Saha-Eggert  relation  if  the  oxygen  enhancement  was  'saturated',  ie  maximised  by 
oxygen  flooding.  They  had  to  reduce  the  02*  primary  current  to  achieve  this  as  too  fast 
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an  erosion  rate  prevented  the  effective  adsorption  of  oxygen. 

Wilson  et  al22  give  the  relative  sputter  rates  for  Si,  GaAs,  InP  and  CMT  in  the  ratios 
1,2,3  and  5.6.  As  the  concentration  of  implanted  oxygen  varies  inversely  as  the  sputter 
rate,  these  figures  are  consistent  with  the  theory  that  in  CMT,  and  possibly  in  InP,  the 
concentration  of  implanted  oxygen  is  insufficient  to  give  optimum  enhancement.  It  may 
be  possible  to  correlate  the  deviation  from  linearity  with  relative  affinities  for  oxygen. 
It  would  be  interesting  to  test  this  theory  by  remeasuring  the  RSFs  with  a  lower  primary 
ion  current  and  oxygen  back-  flooding  though  the  results  would  be  of  more  academic  than 
practical  interest  as  the  resulting  RSFs  would  only  be  applicable  under  such 
experimental  conditions.  If  successful,  however,  it  would  lend  weight  to  the  above 
model  and  encourage  further  theoretical  work  on  the  effects  of  partial  oxygenation 
which  could  result  in  a  practical  formula  for  RSFs. 

CONCLUSIONS 

The  general  agreement,  within  a  factor  of  2,  between  the  two  groups  of  workers, 
despite  widely  differing  techniques,  shows  that  the  scatter  of  points  on  the  log  (RSF) 
versus  E  plot  for  CMT  is  attributable  not  to  the  experimental  error  but  to  more 
fundamental  aspects  of  the  ionising  process.  Furthermore,  the  results  do  not  fit  the 
Saha-Eggert  equation  which  works  well  for  many  other  matrices,  such  as  metals,  glasses 
and  ceramics. 

Despite  decades  of  basic  studies  there  is  still  no  universally  accepted  model  of  the 
ionisation  process  in  sputtering  and  it  is  remarkable  that  although  the  LTC  model  is 
clearly  unsound,  it  has  given  such  good  quantitative  predictions.  The  cascade  mode!  of 
sputtering  is  well  established.  With  simplifying  assumptions,  eg  ignoring  differences  in 
binding  energy,  one  can  predict  ionisation  probabilities  of  similar  form  to  the  Saha- 
Eggert  equation  just  as  thermodynamic  functions  can  be  derived  from  Boltzman 
statistical  mechanics. 

The  comparison  of  the  semi-log  plots  for  the  semiconductors  Si,  GaAs,  InP  and  CMT 
appears  to  show  a  progressive  increase  of  spread  from  the  trend  lines,  excluding  the 
points  for  Mg  and  Be  which  appear  to  be  in  a  class  of  their  own  and  may  well  involve 
a  different  mechanism.  The  sputter  rates  for  these  matrices  increase  in  the  same  order 
and,  since  the  concentration  of  implanted  oxygen  varies  inversely  with  the  sputter  rate, 
it  suggests  that  the  increasing  spread  may  be  associated  with  lower  concentrations  of 
ion-enhancing  oxygen.  If  it  proves  possible  to  correlate  these  deviations  from  linearity 
with  the  affinity  of  the  trace  elements  for  oxygen,  the  prospect  arises  of  deriving  a 
semi-empirical  relationship  which  would  allow  the  quantitative  prediction  of  RSFs  for 
elements  not  yet  measured  as  well  as  resolving  the  few  discrepancies  between  the 
results  of  the  two  laboratories. 
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TABLE  1  COMPARISON  OF  MEASURED  RSFs 


Element 

RSF: 

RSRE 

Wilson 

Lapides 

Log(Wilson) 

minus 

Log  (RSRE) 

Ag 

1.2e21 

1.2e21 

0.00 

A1 

2.0el8 

4.(*el8 

4.0el8 

0.30 

B 

6.8e20 

9.8e20 

5.6e20 

0.16 

Ca 

4.5el8 

3.3e20 

-0.13 

Co 

2.7el9 

4.9el9 

0.26 

Cr 

2.1el9 

1.6el9 

-0.12 

Fe 

4.6el9 

5.9el9 

0.11 

Ga 

1.0el9 

8  0el8 

5.4el8 

-0.10 

In 

1.9el8 

1.8el8 

1.3el8 

-0.02 

Li 

1.3el8 

8.0el7 

1.3el8 

-0.21 

Mg 

7.7el8 

1.3el9 

1.4el9 

0.23 

Mn 

1.8el9 

3.6el9 

0.30 

Na 

8.6el7 

4.4el7 

4.0el7 

-0.29 

Ni 

2.7e20 

2.2e20 

-0.09 

Si 

5.8e20 

8.0e20 

5.4e20 

0.14 

Sn 

3.8e20 

4.0e20 

0.02 

As 

3.7e21 

(1.2e22) 

2.0e22 

0.73 

Cu 

6.0el9 

2.0e21 

3.6e20 

1.52 

Sb 

2.7e20 

1.5e22 

1.75 

|  Au 

(1.5e21) 

|  Pt 

1.8e22 

1  86 

3.6e20 

3.6e20 

1  Ti 

2.7el9 

Zn 


1.5e22 


TABLE  2  FACTORS  USED  IN  CORRELATION  OF  RSFs 


Element 

Ionisation 

Potential 

bVb°(5500K) 

Ag 

7.58 

0.50 

1  A1 

5.99 

0.17 

B 

8.30 

0.17 

Ca 

6.11 

1.82 

Co 

7.86 

0.83 

Cr 

6.77 

0.66 

Fe 

7.87 

1.51 

Ga 

6.00 

0.19 

In 

5.79 

0.23 

Li 

5.39 

0.47 

Mg 

7.65 

1.95 

Mn 

7.44 

0.87 

mam 

5.14 

0.48 

Ni 

7.64 

0.38 

Si 

8.15 

0.72 

Sn 

7.34 

0.60 

As 

9.81 

1.35 

Cu 

7.73 

0.43 

Sb 

8.64 

0.74 

Au 

9.23 

0.47 

Pt 

... 

9.00 

0.51 

Detector 
Factor  f 


Mean 
log  RSF 


21.08 


18.45 


20.91 


18.59 


19.56 


19.26 


19.72 


18.95 


18.27 


18.01 


19.00 


19.41 


17.79 


20.39 


20.83 


20.59 


22.18 


22.26 


20 


40 


60 


Energy  (eV) 


60  100 


Energy  Distributions  Of  Elements  In  CMT 


Log  (RSF) 


22 


pt+  xAu 


x  As 


oAs 


21-1  +A9 


Ionisation  Potential  (E) 


Fig  2  Log  (RSF)  vs  E-RSRE  data  for  CMT 


sv+ 


(dSH)  Boi 


(referred  to  Si+)  (referred  to  As+) 


APPENDIX  A 


Though  the  excitation  energy  imported  by  cascade  impacts  may  vary  widely, 
Boltzmann's  statistical  approach  may  be  used  to  predict  the  most  likely 
population  distribution  between  the  quantum  levels  for  each  type  of  sputtered 
atom.  The  principal  assumption  is  that  despite  the  scatter  of  excitation 
energy  the  average  energy  is  constant. 


Consider  a  number  of  sputtered  atoms  N.  Let  the  number  excited  to  energy 
levels  e1(  e2  ...  effl  be  n, ,  n2  ...  nm. 

The  average  energy  EA  is  given  by 


NEA  =  n1  e1  *  n2  e2  •  •  *  *  ni  ei  *  •  •  +  nm  em 


(1  ) 


also  N  =  n1  +  n2  +  . . .  ♦  ni  .  . .  ♦  nm 


(2) 


Let  w  the  number  of  ways  in  which  the  N  atoms  may  be  distributed  among  the 
sets  n1#  n2  etc  to  nm.  The  number  of  possible  combinations  is: 


ni  •  •  n2\ - n,!  .  .  .nm! 

and  inW**in(N!)  -  £  ^ n { n , !  ) 

Boltzmann's  assumption  is  that  the  most  likely  distribution,  n1 ,  n2,  n3  etc  is 
that  for  which  W  is  a  maximum,  ie  6  in  W  =  0  for  all  small  changes  6ni  in  the 
variables  n1 . 


Therefore  £  6(€nn , !  )  -  6 <Ln  (Nl )  »  0  (3) 

since  N  is  independent  of  such  changes. 

In  Appendix  B  it  is  shown  that  the  expression 
n1  in  n1  -  n^ 

is  a  good  approximation  for  €n  (n^l)  especially  for  larger  values  of  ni . 


A- 1 


Therefore  6  in  (n,))  =  £n  n,  .  fin, 


and  6 (in  n7!  )  -  finn,  .fin,  *  Cnn2 .  &nz  *  ■  •  •  *  ^nnm  •  &nm 


(4) 


In  general,  if  all  the  6ni  terms  in  an  equation  of  the  type 


£  f  (n,)  .  6m,  -  0 


are  independent  then  the  solution  must  be  f(n7)  =  0  for  each  n . .  In  the 
present  case  however,  the  variables  n,  are  also  subject  to  equations  (1)  and 
(2).  Consequently  there  are  (m-2)  independent  variables  and  the  relations. 


6n1  ♦  6n2  ♦  6n3  +  . . .  +  6nm  =  0 


(5) 


and  e^n,  ♦  e2  6n?  +  e3  6n3  +  ...  ♦  em  6nm  =0  (6) 

Multiplying  (5)  by  a  factor  a  and  (b)  by  (3  and  adding  both  to  (4)  gives 

1  tn 

£  (€nn,  +  a  +  (3e,)  .  6n,  =  0  (7) 

f«i 

by  choosing  a  and  P  to  make 
in  n,  ♦  a  ♦  (Je,  =  0 
and  €n  n2  ♦  a  ♦  0e2  =  0 

equation  (7)  is  reduced  to  the  summation  from  i  =  3  to  i  -  m,  ie  (m-2) 
independent  variables.  Consequently  the  generally  solution  is 


A-2 


£n  ni  +  a  +  3e^  =  0 


for  each  of  the  sets  . 


Therefore  ni  =  exp  (-a  -{5  e^  which  can  be  written  as 


n1  =  A  exp  (-ge^ 


This  is  the  standard  derivation  of  the  Maxwell-Boltzmann  distribution  law,  but 
with  the  assumption  of  thermal  equilibrium  replaced  by  the  idea  that  the 
average  excitation  energy  EA  is  independent  of  N,  the  size  of  the  sample. 
While  it  may  accurately  reflect  the  distribution  of  population  between  the 
energy  levels  of  each  individual  sputtered  element,  the  value  of  EA  and 
consequently  3,  and  possible  a,  may  vary  between  elements  because  of  the 
differences  in  finding  energy.  Similarly  account  must  be  taken  of  the  kinetic 
energy  imported  to  the  particles  by  the  cascade  collisions  but  Fig  1  suggests 
that  the  kinetic  energy  distributions  for  different  elements  tend  to  be 
similar . 

Let  the  energy  levels  of  the  neutral  atoms  be  ev  e2  ...  etc  and  those  of  the 
ionised  atoms,  measured  from  the  ground  state  of  the  ion,  be  e\,  e'z  ...  etc  . 

The  partition  functions  B°  and  B',  listed  in  Table  2,  are  defined  as 


B°  «  £  exp  (  -  (3e,)  and  (3'  )  «  £  exp  (  -3  e’.) 
i  i 

the  summations  being  taken  over  the  energy  levels  of  the  atom  and  ion 
respectively. 

Then  the  ratio  of  ions  to  atoms  is: 

£  exp  [  -3(r  ♦  e\)  ]  /  Y,  exP  < 

-  &  /B°  exp  (  -3D 

A-3 
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Sterling's  approximation  to  factorial  n  is 
n!  -  <j2nn  ( n/e)n 

The  following  table  shows  how  good  this  approximation  is  even  for  small 
numbers : 


n 

Factorial  n 

Sterling 

3 

6 

5.8 

5 

120 

118.0 

1 

5040 

4980.4 

10 

3 . 63e6 

3 . 60e6 

25 

1 . 55e25 

1 . 55e25 

100 

9 . 33e1 57 

9 . 32e1 57 

Taking  logarithms,  ln(n! )  =  $ln  (2nn)  +  n  ln(n)  -  n 

As  n  increases  the  first  term  becomes  less  significant,  so  that: 

fin  (n!  )  =*  nfin(n)-n 
as  shown  in  the  following  table: 


n 

fin (n! ) 

nfin  n-n 

10 

15.10 

13.03 

25 

58.00 

55.47 

1  50 

149.5 

145.6 

I  100 

363.7 

f  300 

1415 

1411__  | 
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